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CHAPTER  1 
INTRODUCTION 


A  common  parameter  for  describing  the  accuracy  of  a  weapon  is  the  circular 
probable  error,  generally  referred  to  as  CEP.  CEP  is  simply  the  bivariate  analog 
of  the  univariate  probable  error  and  measures  the  radius  of  a  mean-centered  circle 
which  includes  50%  of  the  bivariate  probability.  In  the  case  of  circular  normal 
errors  where  the  error  variances  are  the  same  in  both  directions,  CEP  can  be 
expressed  as  a  function  of  the  common  miss  distance  standard  deviation.  Also,  CEP 
estimators  based  on  observed  miss  distances  are  easily  formulated  and  can  be  used 
to  construct  confidence  intervals  for  CEP.  In  the  case  of  elliptical  normal 
errors,  CEP  cannot  be  expressed  explicitly  as  a  function  of  the  miss  distance 
standard  deviations.  Here,  one  must  obtain  CEP  by  numerical  methods  or  by  re¬ 
ferring  to  tabular  values.  This  has  led  to  the  development  of  a  number  of 
approximations  by  which  CEP  can  be  expressed  as  a  function  of  the  miss  distance 
standard  deviations.  While  CEP  estimators  based  on  observed  miss  distances  are 
easily  formulated  from  these  approximations,  their  probability  distributions  are 
too  complicated  to  be  useful  for  CEP  confidence  intervals.  In  this  report,  these 
probability  distributions  are  approximated  with  distributions  which  are  more 
practical  for  the  formulation  and  application  of  CEP  confidence  intervals.  Ap¬ 
proximate  CEP  confidence  intervals  are  then  formulated  and  their  accuracy  deter¬ 
mined  through  Monte  Carlo  sampling. 

The  first  part  of  this  report  is  tutorial  in  the  development  of  CEP  and 
discusses  the  commonly  used  approximations  for  the  elliptical  case.  The  develop¬ 
ment  of  approximate  confidence  intervals  begins  with  Chapter  4. 
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CHAPTER  2 

REVIEW  OF  CIRCULAR  CASE 


In  general,  it  will  be  assumed  that  the  errors  in  the  X  and  Y  directions 

are  independent  with  mean  zero  and  variances  a2  and  a2,  respectively.  Under  the 

x  y 

circular  normal  assumption,  o2  =  o2  =  o2  and  the  bivariate  distribution  of  errors 
is  given  by  ” 


fc(x,y)  = 


1  -(xJ 


2nd2 


:)/2a2 


-  00  < 


x,  y 


<  00 


(2.1) 


where  the  subscript  c  denotes  circular.  The  distribution  of  the  radial  miss 
distance  is  derived  by  first  obtaining  the  distribution  of  the  polar  variables 
(R,  0)  where 

X  =  R  cos  0 

Y  =  R  sin  0. 

This  is  found  to  be 

y  _  y  2  /  p  / x2 

8c(r’8)  =  2nd2  e  ,  0  <  r  <  oo,  0  <  0  <  2 n.  (2.2) 

The  distribution  of  R  =  (X2  +  Y2)^  is  now  obtained  from  (2.2)  by  using  the  mar- 
ginal  rule ;  i.e. , 

8c(r)  =  f  8c(r,e)dfl  =  §*  e-  r2/2a2  ,  r  >  0.  (2.3) 

c  c  a 

This  is  the  well-known  Rayleigh  distribution  (see  Lindgren  (1968))  with  cumulative 
distribution  function 

P(R  <  r)  =  Gc(r)  =  1  -  e  "r2/2°2.  (2.4) 

By  definition,  G  (CEP)  =  .5  and  the  solution  of  (2.4)  yield  the  frequently  used 
expression 

CEP  =  [-2  £n( .50) a  =  1.1774  a.  (2.5) 
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Consider  now  that  a  (and  hence,  CEP)  is  unknown  and  is  estimated  from  n 
observed  miss  distances.  These  miss  distances  will  be  designated  as  (x.,y.), 

i  =  1,  .  .  .  ,  n  in  the  X  and  Y  directions,  respectively.  Moranda  (1959)  has  shown 

that  the  maximum  likelihood  estimator  for  a  is 


a 


n  2  2 

X  (*■  +  Y • ) /2n 
i=l 


(2.6) 


and  the  corresponding  estimator  for  CEP  is  simply  CEP  =  1.1774  cr.  This  estimator 

is  biased  for  CEP,  i.e.,  the  expectation  of  CEP  is  not  equal  to  CEP.  However, 
the  bias  is  small  and  the  unbiasing  factor  is  cumbersome.  Therefore,  it  will 
be  retained  in  its  slightly  biased  form. 


To  place  confidence  limits  on  CEP,  it  will  be  necessary  to  examine  the 
probability  distribution  of  CEP.  It  is  well-known  under  normal  theory  (see  Mood 
and  Graybill  (1963))  that  if  { ,  i  =  1,  .  .  .  ,  n  is  a  random  sample  from  a  normal 
population  with  mean  p  and  variance  a2,  then  X(W^  -  p)2/a2  has  a  chi-square 
distribution  with  n  degrees  of  freedom.  One  can  consider  {x^},  i  =  1,  n 

and  {yi}  ,  i  =  1,  .  .  .  ,  n  to  be  a  random  sample  of  size  2n  from  a  normal  popula¬ 

tion  with  mean  zero  and  variance  a2.  Therefore, 


n 


X 

i=l 


2na2 

a2 


2nCEP2 

CEP2 


(2.7) 


2 

where  designates  ,fis  distributed  asH  and  X2n  designates  a  chi-square  proba¬ 
bility  distribution  with  2n  degrees  of  freedom.  The  100  (1  -  a)%  confidence 
limits  are  now  easily  constructed  using  the  probability  statement 


Prob 


2  2nCEP2 

X2n,a/2  CEP^  X2n,l-a/2 


1  -  a 


(2.8) 


In  this  expression,  x^  a  designates  the  100a  percentage  point  for  a  chi-square 

with  u  degrees  of  freedom.  Tabular  values  for  integral  u  can  be  found  in  the 
back  of  most  statistics  texts.  A  more  complete  table  is  found  in  Hald  (1952). 
Manipulating  the  inequality  in  (2.8)  leads  to  the  following  100  (1  -  a)%  con¬ 

fidence  limits  for  CEP: 


A 

CEP 

A 

CEP 

.(x2n,l-<,/2/2n^  ’ 

^X2n,a/2^2n  ^ 

(2.9) 


The  interpretation  here  is  that  one  is  100  (1  -  a)%  confident  that  the  interval 
in  (2.9)  contains  the  population  CEP.  This  formula  is  valid  only  for  the  case 

where  the  errors  are  known  to  be  circular,  i.e.,  the  case  where  a2  =  a2  =  a2. 

>  >  x  y 
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Before  leaving  this  review  of  the  circular  case,  it  will  be  instructive  to 
work  through  an  example.  Suppose  confidence  limits  on  CEP  are  desired  from  the 

ten  round  sample  shown  in  Table  2-1.  We  first  need  to  compute  cr  in  (2.6).  One 
notes  that  the  sum  under  the  radical  in  (2.6)  can  be  expressed  as 


Ex? 

i 

n 


TABLE  2-1.  10  HYPOTHETICAL  MISS  DISTANCES  (FEET) 


42 

-123 

-13 

-12 

-50 

14 

-70 

169 

-191 

-58 

117 

-79 

158 

99 

16 

-18 

101 

170 

27 

65 

The  two  components  are  independent  estimates  of  the 
differ  significantly,  it  would  cast  doubt  on  the 

These  components  will  be  referred  to  as  s2  and  s2 

x  y 


a  = 


'H  +  syV2] 


common  variance  a2.  If  they 
circular  normal  assumption. 

so  that  a  in  (2.6)  becomes 


For  this  example,  one  finds 


s2  =  Ix2/n  =  9565.3 
x  i 

s2  =  Zy2/n  =  9688.5 
y  Ji 

a  =  [(9565.3  +  9688. 5)/2] ^  =  98.12 


CEP  =  1.1774  a  =  115.53. 


To  form  confidence  limits,  the  computations  in  (2.9)  are  required.  Let  us 
consider  95%  limits  so  a.  =  .05  and  the  tabular  values  required  are 

X20 , . 025  =  9,59 
X20 , . 975  =  34-20, 

These  would  both  be  divided  by  2n  =  20  to  form  the  terms  under  the  radical  in 
(2.9).  One  could  also  use  a  table  of  chi-square  percentage  points  divided  by  the 
degrees  of  freedom  here  to  avoid  the  latter  step.  Such  a  table  is  in  Hald  (1952) 
and  provides 
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*20,. 025^°  •Zt796 
X20,.975/20  =  1*7085 


The  95%  confidence  limits  on  CEP  can  now  be  completed  and  are  found  to  be 

/ — 115.53^  ^  115. 53^ \  _  (88  39>  166.82) 

\(1.7085)"2  (.4796) V 


The  units  are  feet,  the  same  as  the  miss  distance  units  in  Table  2-1.  The 
interpretation  is  that  one  is  95%  confident  that  the  true  (or  population)  CEP 
lies  in  the  interval  (88.39,  166.82).  The  result  is  valid  only  if  the  probability 
distribution  of  miss  distances  follows  a  circular  normal  distribution.  Applica¬ 
tion  of  (2.9)  when  the  probability  distribution  is  elliptical  can  lead  to  serious 
errors.  A  discussion  of  the  elliptical  case  begins  with  Chapter  3. 


2-4 


NSWC  TR  83-205 


CHAPTER  3 

CEP  DERIVATION  AND  APPROXIMATIONS  FOR  ELLIPTICAL  ERRORS 


In  the  elliptical  case,  the  error  variances  are  unequal  and  the  bivariate 
distribution  of  errors  is  given  by 


f£(x,y) 


1  -^[(x/CTx)2  +  Cy/CT  )2] 

2nd  ct  e 
x  y 


oo  <  x,y  <  oo 


where  the  subscript  E  denotes  elliptical.  The  distribution  of  the  radial  error 
R  for  this  case  was  derived  by  Chew  and  Boyce  (1961).  They  proceeded,  as  in 
the  circular  case,  by  first  obtaining  the  distribution  of  the  polar  variables 
(R,0).  This  was  found  to  be 


where 


8E(r,6) 


r 

2ttct  o  e 
x  y 


ar2  e(-br2cos20) 


0  <  r  <  oo 

0  <  0  <  2n 


a  = 


oz  +  oz 
-  JL 


(2 ct  a 

v  x  y  / 


b  = 


oz 
-  y 


(2  a  a  V 
'  x  y/ 


(3.1) 


Using  the  marginal  rule,  the  distribution  of  R  was  obtained  by  integrating  g„(r,0) 

Hi 

in  (3.1)  with  respect  to  8  between  0  and  2n.  This  integration  cannot  be  expressed 
in  tractable  form,  so  they  expressed  their  result  in  terms  of  a  modified  Bessel 
function  as 


2 

ggU)  =  o~o~  e  '  3r  Io(br2),  0  <  r  <  oo  .  (3. 2) 

x  y 

In  this  expression,  the  subscript  E  denotes  elliptical  and  I0  is  a  modified  Bessel 
function  of  the  first  kind  and  zero  order,  i,e., 

T  /  \  1  r1  -xcos0 

I0(x)  =  -  /  e  d8. 

7T  ^0 

The  cumulative  distribution  function  for  R  is  denoted  by 

P(R  <  r)  =  GE(r)  =  g£(t)  dt.  (3.3) 
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However,  Gg(r)  cannot  be  expressed  in  tractable  form  because  g^(t)  cannot  be  so 

expressed.  This  means  that  the  radius  of  the  50%  circle  for  the  elliptical  case 
cannot  be  expressed  by  a  simple  formula  as  it  was  in  the  circular  case.  One  has 
to  solve  Gg(CEP)  =  .5  by  numerical  methods  or  by  referring  to  tables  prepared  by 

Harter  (1960),  DiDonato  and  Jarnagin  (1962),  and  others.  To  avoid  using  these 
tables  or  numerical  procedures  for  CEP  evaluation,  a  number  of  approximations  have 
been  developed  over  the  years.  Five  of  the  most  common  are  shown  below;  they  are 
designated  as  CEPX  through  CEP5: 


CEPX 


1.1774 


ct2  +  a2  \  ^ 

-V) 


CEPc 


1.1774 


a  +  a 

je _ y 


cep3 


2  X 


u,  .50 


u  = 


a5  +  a4 

x  y 


a2  +  a2 
JE _ 1 


2 


^2 


CEP4 


cep5 


565 

CT  + 

max 

.612 

a  .  , 
mm  ’ 

0  .  /ct 
min  max 

> 

.25 

667 

CT  + 

max 

.206 

a  .  , 
min  ’ 

CT  .  /ct 
mm  max 

< 

.25 

CEPX  and  CEP2  were  taken  from  Groves  (1961);  CEP3  was  formulated  using  the  chi- 

square  approximation  for  calculating  hit  probabilities  provided  by  Grubbs  (1964). 
It  was  also  derived  independent  of  the  Grubbs  approximation  by  Thomas  and  Taub 
(1978).  CEP4  is  a  piece-wise  linear  combination  of  standard  deviations  which  is 

commonly  used  in  the  missile  community;  CEP5  was  formulated  by  Terzian  (1974) 

using  the  Wilson-Hilferty  approximation  for  calculating  hit  probabilities  provided 
by  Grubbs  (1964).  Plots  of  each  approximation  versus  the  true  CEP  as  a  function 
of  ^min^max  are  s^own  i-n  Figures  3-1  through  3-5.  These  give  a  fairly  good 

indication  of  how  well  each  performs.  It  is  seen  that  CEPX  deteriorates  rapidly 

as  we  depart  from  the  circular  case  (for  which  CEPX  degenerates  to  1.1774  ct)  , 

CEPo  is  reasonably  good  if  the  ratio  ct  /ct  is  not  less  than  about  .2;  CEP3 

appears  good  for  all  ratios,  and  CEP4  and  CEP5  appear  good  to  a  lesser  extent 

for  all  ratios. 
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FIGURE  3-1.  CEPj  APPROXIMATION  VERSUS  TRUE  CEP 


FIGURE  3-2.  CEP2  APPROXIMATION  VERSUS  TRUE  CEP 
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FIGURE  3-3.  CEP3  APPROXIMATION  VERSUS  TRUE  CEP 


CEP,  = 


•565  <Jmax  +-612amin  - 


^min 


a. 


max 


^  .25 


FIGURE  3-4.  CEP4  APPROXIMATION  VERSUS  TRUE  CEP 
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FIGURE  3-5.  CEP5  APPROXIMATION  VERSUS  TRUE  CEP 


To  recap  the  elliptical  case  thus  far,  for  known  values  of  a  and  a  ,  one 

x  y  ’ 

can  compute  the  exact  CEP  by  solving  Gg(CEP)  =  .50  or  compute  an  approximate  CEP 

by  using  an  approximation  such  as  CEPX  to  CEP5.  Consider  now  that  a  and  o  (and 

x  y 

hence  CEP)  are  unknown  and  are  to  be  estimated  from  n  observed  miss  distances. 
These  miss  distances  will  be  designated  as  before,  as  (x^,y^),  i  =  1,  .  .  . ,  n  in 
the  X  and  Y  directions,  respectively.  The  maximum  likelihood  estimators  for  a 
and  0^  pose  no  problem.  They  are  given  by 

In  |  h  -  n 

£xi/nf  ,  sy  =  Oy  ={  £  ypn  (3.4) 

as  shown  in  Lindgren  (1968).  These  estimators  are  slightly  biased  and  as  before, 
they  will  not  be  corrected  due  to  the  cumbersome  nature  of  the  correction  factor. 
They  can  now  be  substituted  for  the  unknown  a’s  in  (3.3)  to  obtain  a  numerical 
estimate  of  CEP  by  solving  Gg(CEP)  =  .50  or  they  can  be  substituted  into  CEP!  to 

CEP5  to  obtain  estimates  of  the  approximate  CEP.  These  latter  estimators  will 

A  A 

be  referred  to  as  CEP^  to  CEP5  and  have  the  appeal  of  being  explicitly  expres¬ 
sible.  Hence,  estimates  of  CEP  for  the  elliptical  case  can  be  rather  easily 
obtained . 
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The  problem  of  obtaining  confidence  limits  for  CEP  in  this  case  is  more 
complex.  The  complexity  is  based  on  the  fact  that  to  find  confidence  limits  for 
a  parameter,  one  needs  information  regarding  the  probability  distribution  of  the 

A 

estimator  for  the  parameter.  As  previously  noted,  CEP  can  be  obtained  by  solving 
Ge(CEP)  =  .50.  Symbol  ically,  we  can  write 

CEP  =  Ge  _1  (.50)  (3.5) 

but  to  obtain  it  requires  recursive  numerical  integration  or  the  use  of  previously 
noted  tables.  Hence,  the  formulation  of  confidence  limits  based  on  this  esti¬ 
mator  holds  little  promise  for  a  practical  solution.  Therefore,  we  shall  consider 
the  formulation  of  confidence  limits  based  on  the  estimators  of  the  approximate 
CEP.  The  distributions  of  these  estimators  are  extremely  complicated  since  they 
involve  radicals  and  linear  combinations  of  sample  variances  and  standard  devi¬ 
ations.  Hence,  these  distributions  were  approximated  and  confidence  limits 
formulated  on  the  basis  of  the  approximate  distributions.  This  development  is 
provided  in  Chapter  4. 
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CHAPTER  A 

APPROXIMATE  DISTRIBUTIONS  OF  CEP  ESTIMATORS 


The  five  estimators  for  CEP  fall  into  two  classes.  One  class  involves  the 
square  root  of  linear  combinations  of  sample  variances  and  the  other  involves 

linear  combinations  of  sample  standard  deviations.  CEP1}  CEP3,  and  CEP5  fall  into 
the  first  class  and  can  be  written  in  the  form 

(s2  +  s2\  ^ 

CEPi  =  Ki  —2~  1  ,  i  =  1,  3,  5  (4.1) 

where 

Kx  =  1.1774 

K3  =(2  Xu,.50/U)h 

Ks  =  [21/3  0  -  2/90)]  3/2  . 

A  A 

CEP2  and  CEP4  fall  into  the  second  class  and  can  be  written  in  the  form 


CEP.  =  ai  s 

1  1  max 

+  a2 

s  .  ,  i 

mm  ’ 

=  2, 

4 

(A. 2) 

where  for 

i  =  2,  ax  = 

a2  = 

1.1774/2 

and  for 

i  =  4,  a1  = 

.565 

and  a2  = 

.612 

when 

s  .  /s  >  .25 

min  max  - 

a  1  = 

.667 

and  a2  = 

.206 

when 

s  .  /s  <  .25. 

min  max 

The  distribution  of  the  square  of  each  estimator  can  be  approximated  by  a 
chi-square  distribution  with  appropriate  degrees  of  freedom.  The  following  is  the 

rationale  for  these  approximations.  The  squares  of  CEPi,  CEP3,  and  CEP5  are 

linear  combinations  of  sample  variances.  Satterthwaite  (1946)  has  shown  that 
one  can  approximate  the  distribution  of  such  linear  combinations  with  a  chi-square 
distribution  with  degrees  of  freedom  chosen  so  the  approximate  distribution  has  a 
variance  equal  to  that  of  the  exact  distribution.  Here,  one  is  approximating  the 
distribution  of  a  linear  combination  of  sample  variances  with  a  chi-square.  A 
natural  extension  is  to  approximate  the  distribution  of  a  linear  combination  of 
sample  standard  deviations  with  a  chi-distribution  (see  Appendix  C) .  Hence, 

the  distributions  of  CEP2  and  CEP4  were  approximated  by  a  chi-distribution  with 

degrees  of  freedom  chosen  so  the  approximate  distribution  has  a  variance  equal 
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to  that  of  the  exact  distribution.  Since  the  square  of  a  chi-variable  has  a 

chi-square  distribution,  the  squares  of  CEP2  and  CEP4  also  have  approximate  chi- 

square  distributions.  The  major  task  is  to  find  the  appropriate  degrees  of  free¬ 
dom  for  each  class  of  estimator. 

First,  one  sees  that  CEPX  has  the  same  form  as  the  maximum  likelihood 
estimator  for  CEP  in  the  circular  case.  In  that  case,  the  degrees  of  freedom 
associated  with  CEP  were  2n.  These  same  degrees  of  freedom  will  be  retained  here 
for  CEPj.  This  will  eventually  show  how  poorly  CEPX  performs  as  an  estimator  for 
CEP  when  the  error  distribution  is  elliptical  vice  circular. 

Next,  consider  CEP3  and  CEP5  of  the  first  class.  Each  has  form 


To  obtain  u',  the  degrees  of  freedom  for  our  chi-square,  we  need  to  equate  the 


variance  of  l)'CEP?/CEP?  with  2o*  (the  variance  of  a  chi-square  with  u’  degrees 
11 

of  freedom)  and  solve  for  o’ .  Now 


CEP?~  2  CEPy 


1  1 


and  the  variance  of  this  expression  is 


i 


this  becomes 


Equating  this  expression  to  2u'  and  solving  for  u*  yields 


(4.3) 


where  u  was  previously  defined  to  be 
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u  = 


k  :  -a* 

a4  +  a4 


(4.4) 


Consider  next,  CEP2  and  CEP4  of  the  second  class.  Each  of  these  estimators 
has  the  form 


CEP.  =  ai  s  +  ao  s 
x  1  max  * 


mm 


i  =  2,  4  . 


To  obtain  0*,  the  degrees  of  freedom  for  the  chi-square  in  this  class,  we  equate 


(u *)'2  CEP. 


the  variance  of 


CEP. 

1 


—  with  0*  -  2 


u*  +  1 


tf) 


(the  variance  of  a  chi-random 


variable  with  U*  degrees  of  freedom)  and  solve  for  u*.  Now, 

(u*)%  CEP.  (u*)^  ( &1  s  +  a2  s  .  ) 

1  _  \  1  max  ^  ram; 


CEP. 

1 


( 


ai  a  +  a o  o 

x  m  o  v  *-  rr 


max  min 

and  the  variance  of  this  expression  is 

U*  faf  V  (  s  \+  a§  V  (  s  .  V\ 

V  _ max/  _  V  mi ny 

/ai  a  +  a2  cr  2 

V  max  z  nun ) 

If  we  denote  the  function  H(x)  as 

+  1 


) 


(4.5) 


Hix)  =l/f 


r 


t!)  ’ 


(4.6) 


then  the  variance  of  a  sample  standard  deviation  based  on  n  observations,  i.e., 
the  variance  of  s^  or  s^  shown  in  (3.4)  can  be  expressed  as 

V(s  )  =  V(s  )  =  CT2  [1  -  H2 (n) ] . 

A  y 

Using  this  notation  in  (4.5)  and  equating  the  latter  to  the  variance  of  a  chi- 
random  variable  yields 


(a?  °x  +  °y)(  1  '  H2(n)j 


(ai  CTX  +  a2  a  ) 


=  1  -  H2  (O*) 


or 


H(u*)  = 


1  -  |  1  -  H2 (n)  | 


(a!  ^  +  gy)^ 

(ai  ax  +  a2  a  ) 


(4.7) 


The  value  of  u*  cannot  be  computed  explicitly  but  is  easily  obtained  using  the 
following  procedure. 
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Evaluate  the  right-hand  side  of  (4.7)  using  estimates  of  a  and  a  obtained 

x  y 

from  n  sample  data  points  and  values  of  aj  and  a2  determined  by  (4.2).  Call 

this  value  w.  Refer  to  Appendix  B  which  contains  tabled  values  of  x  and  H(x)  . 
Enter  the  table  and  find  the  value  of  x  for  which  H(x)  =  w.  This  value  of  x 
is  IK.  An  example  which  incorporates  this  procedure  begins  in  Chapter  5. 

Clearly,  o'  and  u*  may  take  on  fractional  (non-integral)  values.  However, 
this  poses  no  problem.  Although  the  question  of  fractional  degrees  of  freedom 
is  rarely  addressed  in  standard  statistics  texts,  an  extensive  table  of  chi- 
square  percentage  points  with  fractional  degrees  of  freedom  has  been  generated 
by  DiDonato  and  Hageman  (1977).  Also,  one  can  obtain  such  percentage  points 
using  the  MDCHI  subroutine  available  in  IMSL  (1982). 

For  each  of  the  five  CEP  estimators,  it  has  been  shown  that  the  distribution 
of 

U.  CEP? 

CEP^  i  -  1 ,  2 ,  ...  5  (4.8) 

l 

can  be  approximated  by  a  chi-square  distribution  with  LK  degrees  of  freedom  where 

LL  is  either  2n,  o1  ,  or  u*  defined  previously.  Since  the  form  of  a  confidence 

interval  for  a  chi-square  random  variable  is  well-known,  construction  of  con¬ 
fidence  intervals  for  CEP,  using  (4.8),  is  straightforward. 


4-4 


NSWC  TR  83-205 


CHAPTER  5 

APPROXIMATE  CEP  CONEIDENCE  INTERVALS 


An  approximate  100  (1  -  a)%  confidence  interval  for  the  true  (population) 
CEP  can  be  constructed  using  the  probability  statement 


Prob  \  X2 


U.  CEP2 


y  X  X  ^  2 

'u.  ,  a/2  CEP?  V,  l-a/2 


=  1  -  a  . 


(5.1) 


The  subscript  i  is  used  to  indicate  the  approximation  on  which  the  estimate  CEP^^ 
and  the  degrees  of  freedom,  ll,  are  based.  Rewriting  (5.1)  in  terms  of  CEP^^  yields 


CEP. 

i 


r  <  CEP.  < 


CEP. 

i 


fXu.,l-c/2/oi)  1  (XS, ,«/2/ui  ) 


(5.2) 


However,  CEP^  represents  an  approximation  to  the  true  CEP  for  any  i.  Therefore, 

(5.2)  may  be  considered  an  approximate  confidence  interval  for  CEP  and  expressed 
as 


CEP. 

i 


T  <  CEP  < 


CEP. 

l 


Ml-,  l-a/2/ui)^  (xS.,a/2/u1)  ‘ 


(5.3) 


In  the  following  example,  approximate  confidence  intervals  will  be  computed 
for  CEP  using  two  CEP  estimators,  CEP3  and  CEP4.  Using  (3.4),  estimates  of  a 
and  a  can  be  computed  for  the  12  sample  miss  distances  given  in  Table  5-1. 


TABLE  5-1.  12  HYPOTHETICAL  MISS  DISTANCES  (FEET) 


* _  _ y 


-163 

-363 

104 

-56 

-47 

224 

-13 

-61 

-84 

-267 

53 

-85 

93 

383 
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TABLE  5-1.  (Cont.) 


* _  _ y 


197 

-147 

-266 

61 

135 

626 

107 

187 

-112 

-11 

For  these  data,  one  finds 
n 

s2  =  E  x2/n  =  17,505.00 
X  i=l  1 


s 


2 

y 


n 

E  Y ?/n  =  72,191.75 
i=l 


As  previously  defined, 


cep3 


2x5,.50/u 


Since  u  assumes  values  between  1  and  2  inclusively,  a  table  of  the  (2xu  ^q/u)^ 

factor  is  readily  constructed  using  chi-square  percentage  points  taken  from 
DiDonato  and  Hageman  (1977).  Table  5-2  is  a  short  table  that  has  been  prepared 
to  facilitate  computation. 


TABLE  5-2. 


MULTIPLYING  FACTORS  OF 


FOR  CEP3 


APPROXIMATION 


(2  Xu,.50/U)^ 


1.0 

.9538 

1.1 

.9928 

1.2 

1.0258 

1.3 

1.0542 

1.4 

1.0789 

1.5 

1.1005 

1.6 

1.1195 

1.7 

1.1365 

1.8 

1.1516 

1.9 

1.1652 

2.0 

1.1774 

_(gx  *  °v) 

a4  +  a4 

x  y 
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The  12  sample  miss  distances  in  Table  5-1  are  used  to  obtain  estimates  of  u 
and  u?  given  by  u  and  u1  below: 


o  =  4  — a —  =  1.46  o !  =  nu  =  17.52 

s^  +  s^ 


s^  +  s  ’ 
x  y 


Interpolation  in  Table  5-2  gives  ^2  50^°)^  = 


1.0919  so  that 


CEP3  =  (1.0919)  (211.77)  =  231.23 


To  form  confidence  limits,  the  computations  in  (5.3)  are  needed.  If  95%  limits 
are  considered,  the  chi-square  tabular  values  for  a  =  .05  are  obtained  via 
interpolation  in  DiDonato  and  Hageman  (1977)  and  are 


7.91 


Xu' , .025 


The  approximate  95%  confidence  interval  for  CEP  based  on  approximation  3  is, 
therefore 


or 


(174.14  <  CEP  <  344.13)  . 


The  interpretation  here  is  that  one  is  approximately  95%  confident  that  the  true 
CEP  lies  in  the  computed  interval. 

Let  us  next  consider  the  approximate  95%  confidence  interval  obtained  by 

using  CEP4.  Before  using  estimator  4,  one  must  compute  the  ratio  c  =  s  / smax 

to  determine  which  half  of  the  piece-wise  approximation  should  be  used.  In  this 
case,  c  =  .49  so  that 


CEP4  =  .565  s 


max 


+  .612  s  .  =  232.78. 

min 


This  is  reasonably  close  to  the  231.23  obtained  for  CEP3. 


To  determine  u*,  an  estimate  of  U*,  evaluate 
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H  (n)  may  be  determined  using  tabled  values  of  the  gamma  function  provided  in 
the  National  Bureau  of  Standards  Applied  Mathematics  Series  document  by  Salzer 
(1951).  An  abbreviated  version  is  given  in  Appendix  A.  H(n)  may  also  be  read 
directly  from  the  table  provided  in  Appendix  B. 

Now 


H2 (n)  = 


.9592 


and 

A 

H(u*)  =  .9888  . 

Entering  Appendix  B  with  .9888  and  interpolating  between  .9887  and  .9890  yields 
u*  =  22.17.  From  DiDonato  and  Hageman  (1977),  obtain 


.975 


37.00 


via  interpolation.  An  approximate  95%  confidence 
fourth  estimator  is,  therefore,  given  by 


232.78 


(37.00/22 .17)^ 


<  CEP  < 


232.78 


(11. 10/22. 17)^ 


interval  for  CEP  using  the 


or 


(180.19  <  CEP  <  328.98). 

One  notes  that  these  two  intervals  are  different.  Had  the  other  three 
estimators  been  used  to  construct  confidence  intervals,  they  too  would  have 
been  different.  We  now  have  the  problem  of  deciding  which  estimator  to  use  for 
constructing  confidence  intervals.  This  will  be  discussed  in  the  next  chapter. 

Before  leaving  this  chapter,  something  needs  to  be  said  regarding  the  as¬ 
sumption  of  zero  means.  Throughout  the  development  in  the  report,  it  has  been 
assumed  that  the  errors  have  zero  mean  in  both  directions.  An  error  has  been 
assumed  to  be  a  miss  distance  from  a  target,  and  zero  mean  implies  there  is  no 
bias,  i.e.,  the  target  coincides  with  the  distribution  mean.  There  are  applica¬ 
tions  where  the  errors  are  not  miss  distances,  per  se,  but  deviations  from  the 
mean  impact  point.  This  occurs  when  there  is  bias  in  either  or  both  directions  or 
when  there  is  no  target,  i.e.,  the  firings  are  conducted  to  estimate  dispersion 
without  regard  to  a  target.  In  either  case,  the  impact  distribution  is  no  longer 
centered  on  the  target  but  on  an  unknown  point  (p^,  p^) ,  and  CEP  is  the  radius 

of  the  50%  circle  which  is  centered  on  this  point  vice  the  target.  To  apply  the 
methodology  in  this  report  to  these  cases,  the  squares  of  s  and  s  in  (3.4)  must 
be  modified  to  read  x  ^ 
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n 


s2  =  a2 
x  x 


E  (xi  -  x): 

i=l 

n  -  1 


E  x?  ■  (E  xi)2/n 

i=l  i=l 


n  -  1 


s2  =  a2 

y  y 


E  (y±  -  y): 

1=1^ _ _ _ 

n  -  1 


n 


E  y-  -  (E  yi)2/n 

i=l  i=l 


n  -  1 


(5.4) 


In  (5.4),  x  and  y  are  the  averages  of  the  impact  point  locations  in  the  x  and  y 
directions,  respectively.  In  addition,  the  degrees  of  freedom  associated  with 
each  estimator  must  be  reduced  by  replacing  n  with  n  -  1  in  each  as  follows: 


ESTIMATOR 
CEPi 
CEP  3  ' 

► 

cep5 

J 

cepE 

CEP<i 


MODIFIED  d.f.  FOR  NON-ZERO  MEANS 


2  (n  -  1) 


(n  -  l)u 


H(u*)  = 


(  o  )  (aiCTx  + 

1  -I  1  -  H2  (n  -  1)J  - - - 


(aiQ  +  aoQ  ) 
1  x  *  y 


This  reduction  in  degrees  of  freedom  means  a  slight  reduction  in  precision.  In 
effect,  it  requires  n  +  1  deviations  from  the  mean  to  provide  the  same  precision 
as  n  miss  distances. 
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CHAPTER  6 

MONTE  CARLO  VERIFICATION 


To  ascertain  if  the  formulated  approximations  produced  confidence  intervals 
with  confidence  close  to  1  -  a,  a  Monte  Carlo  simulation  was  written  for  the 
CDC  6700.  This  would  also  provide  a  means  of  comparing  the  estimators  with 
respect  to  their  confidence  level  and  their  expected  confidence  interval  length. 
One  replicate  of  the  simulation  entailed  generating  a  sample  of  n  miss  distances 

from  a  bivariate  normal  distribution  with  zero  mean  and  variances  o ^  and  cr^ 

x  y ' 

These  n  values  were  then  used  to  compute  a  CEP  confidence  interval  using  the 

five  estimators  CEPi  to  CEP5.  The  length  of  each  interval  was  computed,  and  for 

each  interval,  it  was  determined  whether  the  interval  contained  the  true  (popula¬ 
tion)  CEP.  This  process  was  replicated  N  =  10,000  times.  The  proportion  of 
replicates  in  which  the  confidence  interval  contained  the  true  CEP  provided 
an  estimate  of  the  confidence  associated  with  each  estimator,  and  the  average 
interval  length  for  each  provided  an  estimate  of  the  expected  length. 

The  values  of  the  parameters  used  in  the  simulation  are: 


n  -  number  of  miss  distances  =  sample  size  =  5,  10,  20 
N  =  number  of  replicates  =  10,000 


1  -  a  -  nominal  confidence  level  =  .95 


c  =  a  .  /ct 

mm  max 


1.0,  .75,  .50,  .35,  .20,  .05. 


It  was  necessary  to  run  a  large  number  of  replicates  to  ensure  results  with 
reasonable  precision.  The  10,000  replicates  used  provided  estimates  of  confidence 
with  an  error  of  less  than  .01  with  probability  .95.  Because  of  the  large  N, 
the  sample  sizes  were  restricted  to  small  values  to  keep  the  computer  time 
within  bounds.  A  nominal  value  of  a  —  .05  ( —  ^  a  nominal  confidence  level  of 

.95)  was  used  in  the  construction  of  all  confidence  intervals  within  the  simula¬ 
tion.  It  follows  that  if  the  confidence  intervals  were  exact  (vice  approximate), 
the  Monte  Carlo  confidence  estimates  should  be  within  .01  of  .95  or  .95  ±  .01 
with  probability  .95.  Hence,  any  confidence  estimate  which  departs  seriously 
from  .95  ±  .01  will  reflect  poorly  on  the  estimator  which  produced  it.  One 

additional  comment  is  required  before  discussing  the  results.  The  simulated 
confidence  levels  depend  only  on  the  ratio  of  the  sigmas  so  that  it  was  not  neces¬ 
sary  to  vary  both  standard  deviations.  The  larger  was  designated  a and  set  equal 

to  unity  while  the  smaller  was  designated  and  set  equal  to  the  ratio  c.  While 

the  average  confidence  interval  lengths  are  dependent  on  both  sigmas,  they  were 
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only  constructed  for  a  =  c,  a  =1 

x  5  y 

purposes.  If  such  lengths  are  needed  for  values  of  a  and  O  other  than  c  and  1, 

x  y  ’ 


This  is  all  that  is  required  for  comparison 

for  values  of  a  and  a  other  than  c  and  1, 
x  y  * 

they  can  be  obtained  by  multiplying  the  tabled  entry  for  appropriate  c  by  o  = 

a  .  y 


max 


The  results  of  the  simulations  are  set  out  in  Tables  6-1  and  6-2.  Let  us 

first  discuss  the  simulated  confidence  levels  shown  in  Table  6-1.  One  first 
notes  that  all  five  estimators  provide  confidence  within  (or  nearly  within)  the 
sampling  variations  (±.01)  of  .95  when  c  =  1.  This  is  the  circular  case  and 

A 

all  the  estimators  except  CEP5  degenerate  properly  to  the  maximum  likelihood 

estimator  for  CEP.  Hence,  the  result  that  all  do  well  for  c  =  1  is  not  unex¬ 
pected.  Next,  one  notes  that  as  c  departs  from  unity  (that  is,  as  the  impact 

distribution  departs  from  circularity) ,  the  confidence  associated  with  CEPi 

departs  seriously  from  .95.  For  example,  with  n  =  10  and  c  =  .20  (5  to  1  ratio 

of  the  sigmas)  the  confidence  associated  with  CEPi  is  only  .689.  This  means  that 
if  one  were  to  use  circular  theory  to  construct  a  95%  confidence  interval  for 

CEP  when  the  distribution  was  elliptical  with  a  5  to  1  ratio  of  the  sigmas,  his 

interval  would  have  confidence  of  less  than  .7!  This  rules  out  CEPi  unless  one 

is  nearly  certain  that  the  impact  distribution  is  circular  normal.  This  result 
is  also  not  unexpected,  but  it  does  quantify  how  poorly  the  circular  estimator 
performs  in  the  elliptical  case. 

In  general,  the  others  do  reasonably  well  unless  c  is  small.  One  notes  this 
especially  for  CEP2  when  c  =  .05;  the  confidence  falls  from  .930  for  n  =  5  to 
.912  for  n  =  20.  It  would  continue  to  decrease  as  n  increases  duq  to  the  error 

A 

in  CEP  approximation  for  small  c  (see  Figure  3-2).  The  distribution  of  CEP2 
becomes  more  concentrated  about  the  approximate  CEP  as  n  increases.  If  the 
approximation  is  in  serious  error  (which  it  is  for  small  c),  then  the  distribution 
is  concentrated  about  the  wrong  value.  The  simulation  was  run  for  n  =  100  at 
c  =  .05  with  a  resulting  confidence  estimate  of  only  .714.  We  see  the  same 

behavior  at  larger  c  values  but  to  a  lesser  extent.  For  example,  at  c  =  .35, 
the  confidence  estimate  is  .942  for  n  =  20  but  dips  to  .922  for  n  =  100.  The 

upshot  here  is  that  CEP2  would  be  a  problem  for  small  c  or  even  moderate  c  if 

A 

the  sample  is  large  enough  With  regard  to  CEP3,  one  sees  that  small  values  of 

c  pose  no  problem.  In  fact,  the  confidence  for  CEP3  is  asymptotic  to  .95  at 

c  —  0 .  Also,  for  values  of  c  around  .5,  the  confidence  estimates  are  slightly 
higher  than  .95.  It  tends  to  peak  out  at  about  .97.  Selected  runs  for  n  =  100 
show  that  this  result  changes  very  little  with  n.  There  is  a  slight  price  to 
pay  for  this  extra  confidence,  and  this  will  be  addressed  when  Table  6-2  is 

discussed.  CEP4  provides  confidence  close  to  .95  for  all  values  of  c  except 

those  where  CEP4  departs  from  CEP  (see  Figure  3-4).  At  those  values,  there  is 

a  reduction  in  confidence  which  increases  with  n  but  not  as  severely  as  for 

CEP2.  The  performance  of  CEP5  is  not  poor  with  respect  to  confidence.  However, 
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TABLE  6-1. 

SIMULATED 

n  =  5 

CONFIDENCE 

LEVELS 

A 

c 

CEPj. 

cep2 

cep3 

cep4 

cep5 

1.0 

.950 

.947 

.963 

.946 

.963 

.75 

.941 

.947 

.965 

.945 

.965 

.50 

.894 

.941 

.968 

.944 

.967 

.35 

.830 

.937 

.963 

.939 

.961 

.20 

.753 

.932 

.952 

.923 

.950 

.05 

.714 

.930 

.950 

.936 

.948 

n  = 

10 

C 

CEPi 

cep2 

cep3 

cep4 

cep5 

1.0 

.947 

.945 

.955 

.944 

.955 

.75 

.935 

.944 

.958 

.943 

.959 

.50 

.876 

.941 

.967 

.943 

.966 

.35 

.789 

.939 

.967 

.941 

.965 

.20 

.689 

.939 

.959 

.931 

.955 

.05 

.640 

.931 

.951 

.943 

.948 

A 

a  = 

20 

C 

CEPi 

cep2 

CEP^ 

cep4 

cep5 

1.0 

.952 

.952 

.956 

.951 

.956 

.75 

.938 

.951 

.960 

.951 

.960 

.50 

.858 

.945 

.970 

.948 

.969 

.35 

.724 

.942 

.970 

.947 

.968 

.20 

.567 

.946 

.961 

.937 

.955 

.05 

.506 

.912 

.950 

.946 

.946 
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TABLE  6-2.  AVERAGE  CONFIDENCE  INTERVAL  LENGTHS 


n  =  5 


A 

A 

c 

CEPt 

CEP2 

CEP3 

CEP4 

CEP, 

1.0 

1.213 

1.160 

1.305 

1.145 

1.317 

.75 

1.071 

1.028 

1 . 175 

1.014 

1.186 

.50 

.950 

.923 

1.118 

.918 

1.131 

.35 

.897 

.889 

1.129 

.917 

1.144 

.20 

.856 

.876 

1 . 149 

.985 

1.168 

.05 

.841 

.912 

1 . 176 

1.054 

1.196 

n  = 

-A 

10 

c 

CEPjl 

cep2 

CEP* 

CEP* 

CEP  5 

1.0 

.792 

.776 

.817 

.768 

.823 

.75 

.698 

.685 

.735 

.676 

.741 

.50 

.622 

.614 

.697 

.601 

.705 

.35 

.587 

.586 

.693 

.586 

.702 

.20 

.564 

.573 

.694 

.636 

.706 

.05 

.553 

.579 

.695 

.662 

.707 

n  = 

20 

A 

C 

CEP* 

cep2 

cep3 

cep4 

CEP* 

1.0 

.537 

.533 

.  545 

.531 

.549 

.75 

.475 

.472 

.492 

.467 

.496 

.50 

.423 

.421 

.465 

.411 

.470 

.35 

.400 

.401 

.460 

.393 

.466 

.20  . 

.384 

.388 

.455 

.429 

.462 

.05 

.378 

.387 

.453 

.441 

.461 
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there  is  little  rationale  for  its  use  in  constructing  confidence  intervals.  The 
Wilson-Hilferty  approximation  avoids  the  use  of  chi-square  percentage  points  for 

fractional  degrees  of  freedom  in  forming  CEP4.  However,  they  are  needed  in  the 

computation  of  the  interval,  so  little  effort  is  saved. 

Let  us  now  discuss  the  average  confidence  interval  lengths  in  Table  6-2.  As 

previously  noted,  these  lengths  depend  on  both  a  and  a  but  were  computed  only 

x  y 

for  values  of  the  ratio  of  cr  to  a  for  comparison  purposes.  However,  we  need 

x  y 

not  compare  all  five  since  some  were  eliminated  as  viable  candidates  in  our 

discussions  of  Table  6-1.  CEP}  was  eliminated  because  its  confidence  eroded 

seriously  as  c  departed  from  unity.  CEP2  had  a  less  serious  but  similar  problem, 

and  CEP5  was  eliminated  because  it  offered  no  improvement  over  CEP3  and  only 

a  slight  reduction  in  computation.  This  leaves  only  CEP3  and  CEP4  to  discuss 

here.  One  notes  that  the  average  lengths  are  uniformly  less  for  CEP4  than  for 

CEP3.  At  mid  values  of  c,  this  is  due  in  part  to  the  inflated  confidence 

inherent  in  the  approximation  of  the  distribution  of  CEP3,  i.e.,  the  higher  the 

confidence,  the  longer  the  confidence  length.  However,  not  all  of  the  difference 
in  length  can  be  attributed  to  higher  confidence.  A  study  by  Taub  and  Thomas 

(1982)  shows  that  the  variance  of  CEP4  is  less  than  the  variance  of  CEP3,  and 

this  is  the  primary  reason  for  the  difference  in  length.  Even  so,  CEP4  suffers 

from  the  bias  caused  by  the  error  in  approximation  shown  in  Figure  3-4.  This  has 

A 

an  effect  on  the  confidence  level  associated  with  CEP4  for  some  values  of  c,  but 
it  would  not  be  appreciable  for  small  n. 

In  summary,  we  can  state  that  the  logical  choice  lies  between  CEP3  and 
CEP4.  The  third  holds  for  all  values  of  c,  regardless  of  n,  and  is  easy  to  imple¬ 
ment.  The  fourth  offers  somewhat  shorter  confidence  lengths  but  it  is  cumbersome 
to  implement  and  would  provide  a  reduced  confidence  level  for  some  values  of  n 
and  c . 

It  would  be  highly  desirable  to  have  a  CEP  estimator  with  a  variance  as 

small  or  smaller  than  CEP4  which  would  have  negligible  bias  for  all  0  <  c  <  1, 

and  which  would  avoid  the  cumbersomeness  of  a  piece-wise  linear  function.  The 
authors  have  several  ideas  along  this  line  and  hope  to  explore  their  merits  in 
the  near  future. 
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APPENDIX  A 

TABLE  OF  T(n)  AND  r(n  +  h)  FOR  n  =  1 ,  ...  100 
Excerpted  from  National  Bureau  of  Standards 
Applied  Mathematics  Series  -  16 
by 


Herbert  E.  Salzer 
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0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


n! 


1.  00000  00000  00000  (0) 

1.  00000  00000  00000  (0) 

2.  00000  00000  00000  (0) 

6.  00000  00000  00000  (0) 

2.  40000  00000  00000  ( 1 ) 

1.  20000  00000  00000  (2) 

7.  20000  00000  00000  (2) 

5.  04000  00000  00000  (3) 

4.  03200  00000  00000  (4) 

3.  62880  00000  00000  (5) 

3.  62880  00000  00000  (6) 

3.  99168  00000  00000  (7) 

4.  79001  60000  00000  (8) 

6.  22702  08000  00000  (9) 

8.  71782  91200  00000  (10) 

1.  30767  43680  00000  (12) 

2.  09227  89888  00000  (13) 

3.  55687  42809  60000  (14) 
6.  40237  37057  28000  (15) 
1  21645  10040  88320  (17) 

2.  43290  20081  76640  (18) 

5.  10909  42171  70944  (19) 

1.  12400  07277  77608  (21) 

2.  58520  16738  88498  (22) 

6.  20448  40173  32394  (23) 

1.  55112  10043  33099  (25) 
4  03291  46112  66056  (26) 

1.  08888  69450  41835  (28) 

3.  04888  34461  17139  (29) 
8.  84176  19937  39702  (30) 

2.  65252  85981  21911  (32) 
8.  22283  86541  77923  (33) 

2.  63130  83693  36935  (35) 
8.  68331  76188  11886  (36) 

2.  95232  79903  96041  (38) 

1.  03331  47966  38614  (40) 

3.  71993  32678  99012  (41) 

1.  37637  53091  22635  (43) 

5.  23022  61746  66011  (44) 

2.  03978  82081  19744  (46) 

8.  15915  28324  78977  (47) 

3.  34525  26613  16381  (49) 

1.  40500  61177  52880  (51) 

6.  04152  63063  37384  (52) 

2.  65827  15747  88449  (54) 

1.  19622  22086  54802  (56) 
5.  50262  21598  12089  (57) 

2.  58623  24151  11682  (59) 
1  24139  15592  53607  (61) 
6.08281  86403  42676  (62) 

3.  04140  93201  71338  (64) 


r  (n+J) 

tl 

n! 

r  (n+*) 

1.  7724 

539 

(0) 

50 

3.  04140 

93201 

71338 

(64) 

4  2904 

629 

(63) 

8.  8622 

693 

(-D 

51 

1.  55111 

87532 

87382 

(66) 

2.  1666 

838 

(65) 

1.  3293 

404 

(0) 

52 

8.  06581 

75170 

94388 

(67) 

1.  1158 

421 

(67) 

3.  3233 

510 

(0) 

53 

4.  27488 

32840 

60026 

(69) 

5.  8581 

712 

(68) 

1.  1631 

728 

(1) 

54 

2.  30843 

69733 

92414 

(71) 

3.  1341 

216 

(70) 

5.  2342 

778 

(1) 

55 

1.  26964 

03353 

65828 

(73) 

1.  7080 

963 

(72) 

2.  8788 

528 

(2) 

56 

7.  10998 

58780 

48635 

(74) 

9.  4799 

344 

(73) 

1.  8712 

543 

(3) 

57 

4  05269 

19504 

87722 

(76) 

5.  3561 

629 

(75) 

1.  4034 

407 

(4) 

58 

2  35056 

13312 

82879 

(78) 

3.  0797 

937 

(77) 

1  1929 

246 

(5) 

59 

1  38683 

11854 

56898 

(80) 

1.  8016 

793 

(79) 

1  1332 

784 

(6) 

60 

8.  32098 

71127 

41390 

(81) 

1.  0719 

992 

(81) 

1  1899 

423 

(7) 

61 

5.  07580 

21387 

72248 

(83) 

6.  4855 

951 

(82) 

1.  3684 

337 

(8) 

62 

3.  14699 

73260 

38794 

(85) 

3.  9886 

410 

(84) 

1.  7105 

421 

(9) 

63 

1.  98260 

83154 

04440 

(87) 

2.  4929 

006 

(86) 

2.  3092 

318 

(10) 

64 

1.  26886 

93218 

58842 

(89) 

1  5829 

919 

(88) 

3.  3483 

861 

(11) 

65 

8.  24765 

05920 

82471 

(90) 

1.  0210 

298 

(90) 

5.  1899 

985 

(12) 

66 

5.  44344 

93907 

74431 

(92) 

6.  6877 

450 

(91) 

8.  5634 

974 

(13) 

67 

3.  6471  1 

10918 

18869 

(94) 

4.  4473 

504 

(93) 

1.  4986 

121 

(15) 

68 

2.  48003 

55424 

36831 

(96) 

3.  0019 

615 

(95) 

2.  7724 

323 

(16) 

69 

1  71122 

45242 

81413 

(98) 

2.  0563 

436 

(97) 

5.  4062 

430 

(17) 

70 

1.  19785 

71669 

96989 

(100) 

1.  4291 

588 

(99) 

1  1082 

798 

(19) 

71 

8.  50478 

58856 

78623 

(101) 

1.  0075 

570 

(101) 

2.  3828 

016 

(20) 

72 

6.  12344 

58376 

88609 

(103) 

7.  2040 

324 

(102) 

5.  3613 

036 

(21) 

73 

4.  47011 

54615 

12684 

(105) 

5.  2229 

235 

( 104) 

1  2599 

063 

(23) 

74 

3.  30788 

54415 

19386 

(107) 

3.  8388 

487 

(106) 

3.  0867 

705 

(24) 

75 

2.  48091 

40811 

39540 

(109) 

2.  8599 

423 

(108) 

7.  8712 

649 

(25) 

76  ; 

1.  88549 

47016 

66050 

(111) 

2.  1592 

564 

(110) 

2.  0858 

852 

(27) 

77 

1.  45183 

09202 

82859 

(113) 

1.  6518 

312 

(112) 

5.  7361 

843 

(28) 

78 

1  13242 

81178 

20630 

(115) 

1.  2801 

692 

(114) 

1.  6348 

125 

(30) 

79 

8.  94618 

21307 

82975 

(116) 

1.  0049 

328 

(116) 

4.  8226 

969 

(31) 

80 

7.  15694 

57046 

26380 

(118) 

7.  9892 

157 

(117) 

1.  4709 

226 

(33) 

81 

5.  79712 

60207 

47368 

(120) 

6.  4313 

187 

(119) 

4  6334 

061 

(34) 

82 

4  75364 

33370 

12842 

(122) 

5.  2415 

247 

(121) 

1  5058 

570 

(36) 

83 

3.  94552 

39G97 

20659 

(124) 

4.  3242 

579 

(123) 

5.  0446 

209 

(37) 

84 

3  31424 

01345 

65353 

(126) 

3.  6107 

553 

(125) 

1.  7403 

942 

(39) 

85 

2.  81710 

41143 

80550 

(128) 

3.  0510 

883 

(127) 

6.  1783 

994 

(40) 

86 

2.  42270 

95383 

67273 

(130) 

2.  6086 

805 

(129) 

2.  2551 

158 

(42) 

87 

2.  10775 

72983 

79528 

(132) 

2.  2565 

086 

(  131) 

8.  4566 

842 

(43) 

88 

1.  85482 

64225 

73984 

(134) 

1.  9744 

450 

(133) 

3.  2558 

234 

(45) 

89 

1.  65079 

55 1 60 

90846 

(136) 

1.  7473 

838 

(135) 

1.  2860 

502 

(47) 

90 

1.  48571 

59644 

81761 

(138) 

1.  5G39 

085 

(137) 

5.  2085 

035 

(48) 

91 

1.  35200 

15276 

78403 

(HO) 

1.  4153 

372 

(139) 

2.  1615 

290 

(50) 

92 

1  24384 

14054 

64131 

(  142) 

1  2950 

336 

( 141  ) 

9.  1864 

981 

(51) 

93 

1  15677 

25070 

81642 

(  144) 

1.  1979 

060 

(  143) 

3.  9961 

267 

(53) 

94 

1.  08736 

6 1 566 

56713 

(146) 

1  1200 

422 

(  145) 

1.  7782 

764 

(55) 

95 

1.  03299 

78488 

23906 

(148) 

1.  0584 

39S 

(  I47i 

8.  0911 

574 

(56) 

96 

9.  91677 

93487 

09497 

(149) 

1.  0108 

100 

v  14<P 

3.  7623 

882 

(58) 

97 

9.  01927 

59682 

48212 

(  151 ) 

9.  7543 

169 

t  150 

1.  7871 

344 

(60) 

98 

9.  42689 

04488 

83218 

(153) 

9.  5104 

590 

(  1 52 ) 

8.  6676 

018 

(61) 

99 

9.  332G2 

15443 

94415 

(155) 

9.  3678 

021 

(  151) 

4.  2904 

629 

(03) 

100 

9.  33262 

15443 

94415 

(157) 

9.  3209 

631 

<  1 56' 
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APPENDIX  B 

TABLES  OF  x  AND  H(x)  FOR  x  =  . 1  TO  400 


8837  3.80  .9370  5.65  .9569 
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APPENDIX  C 

THE  CHI  DISTRIBUTION 
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A  discussion  of  the  chi-square  probability  distribution  can  be  found  in  most 
probability  and  statistics  textbooks.  It  is  a  special  case  of  the  more  general 
gamma  distribution  and  its  density  function  is  given  by 


£XM 


1  xn/2  "  1  -x/2 

r(§)  2  n/2 


0  <  X  <  0° 


Rarely,  however,  is  the  CHI  distribution  discussed  in  detail.  Because  this  dis- 
tribution  is  critical  to  the  development  of  the  approximate  distributions  for 

CEP2  and  CEP4,  the  following  is  a  brief  introduction  to  the  CHI  distribution. 

If  X  is  defined  as  a  chi-square  random  variable,  then  Y  =Vx  is  distributed 
as  a  CHI  random  variable  whose  density  function  is  written  as 


fy(y)  =  - 

r(S)2 

The  mean  of  Y  is  given  by 


n  -  1 


-x_ 

2 


-  1 


0  <  x  <  co 


E(Y)  = 


v( 


n  +  1 
2 


% 


and  the  variance  of 

V (Y)  =  n 

In  general,  the  rt^L 


Y  is  given  by 


moment  of  Y  is 


given  by 


E(Yr) 
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If  s  is  the  sample  standard  deviation  with  u  degrees  of  freedom,  then  the  previous 
results  can  be  used  to  establish  that 


These  results  are  needed  to  determine  the  degrees  of  freedom  associated  with 
CEP2  and  CEP4.  Additional  information  on  the  CHI  distribution  may  be  found  in 
Krutchkoff  (1970). 
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